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THERMODYNAMICS OF QCD AT HIGH TEMPERATURE a 
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Department of Physics, The Ohio State University, Columbus, OH 43210, USA 

A hierarchy of effective field theories is used to separate the contributions from 
different momentum scales and to calculate the free energy of QCD at high tem- 
perature in powers of the coupling constant up to order g 5 . The behavior of the 
perturbative series will also be discussed. 



The partition function of QCD at finite temperature is obtained from the 
Lagrangian £q C d = (1/4)G£„G£„ + q-y^D^q, with gauge coupling constant g. 
We consider the effective theory that results from integrating out the nonstatic 
modes. Such a theory, which is called Electrostatic QCD (EQCD), is made out 
of only the static bosonic modes and 3 dimensional. The free energy density 
is Fqcd = T(fs - logZEQCD/V^), where Zeqcd is the partition function of 
EQCD obtained from the Lagrangian 

£eqcd = \G%G1 + \{D t A a ) 2 + \m\A% + ~\A 4 + S£ EQCD (1) 

whith effective gauge coupling constant <?£. <5£eqcd represents an infinite 
series of non-renormalizable terms. The effects of the fermions are incorporated 
into the effective parameters. By power counting, the effective parameters at 
leading order are: fs ~ T 3 , m% ~ g 2 T 2 , g 2 E ~ g 2 T, and Xe ~ g 4 T. 

The magnetostatic fields remain massless; therefore, we can go further 
in separating the different scales of QCD at high temperature by integrating 
out Aq. We obtain an effective theory of EQCD which is called magnetostatic 
QCD (MQCD). The free energy density of QCD can be written as F = T(f E + 
/m + fa), where f G — - logZ M QCD/^, and -Zmqcd is the partition function 
of MQCD obtained from the Lagrangian £mqcd = (l/4)(r?.Ggf + ^-Cmqcd • 
The gauge coupling constant is gM and <5£mqcd represents an infinite series 
of non-renormalizable terms. 

Again, we can use power counting to identify the order of the leading con- 
tribution to the MQCD parameters: ju ~ "m 3 E ~ g 3 T 3 and g 2 M ~ gs ~ g 2 T. 
\e contributes to the free energy only at order g 6 ; so, if we are interested in the 
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free energy at lower order, we can ignore A e ■ Similarly the non-renormalizable 
terms of EQCD can also be omitted. 

The leading contribution to fa can be obtained by realizing that the only 
parameter with dimensions involved in MQCD is qm and therefore the leading 
contribution is of order 

9 m ~ 9 6 T 3 - Since we are interested in computing the 
free energy of QCD up to order g 5 we can ignore the contribution from fa', it 
will be enough to compute Je and /a/. We conclude that the free energy of 
QCD up to order g 5 is given by F = T[f E (T, g; A E ) + fM(m%, g E ) A E )], where 
As is a factorization scale that separates the scales T and gT. Therefore, we 
have to determine f E , m E , g E , and fu ■ Calculations for SU (A^Vgauge theory 
with rif fermions and results in analytical form are detailed inH; here, I will 
just state the results in numerical form for SU(3) with 3 fermions. 

The effective mass tue is computed by matching the electrostatic screening 
mass for QCD and EQCD. 



m| 



0.93 T 2 

47T 



1 



0.26+ 1.43 log^)£ 



(2) 



At this order in g 2 , there is no dependence on the factorization scale Ag. 

The parameter Je is determined by calculating the free energy in both full 
QCD and EQCD, and matching the two results. 



1 - 0.90 



<£_ 

47T 



18.3 + 6.91 log ^ - 1.30 log 
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2ttT 



<£_ 

-in 



(3) 



where g is the coupling constant in the MS renormalization scheme at the 
scale /x. We have used the renormalization group equation of the coupling 
constant to shift the scale of the running coupling constant to an arbitrary 
renormalization scale /i. 

Through order g 5 , ju is proportional to the logarithm of the partition 
function for EQCD: Jm = — logZEQCD/^- Now, we have to consider the 
cotnribution to log -Zeqcd of orders g 3 , g , and g b which are given by the sum 
of 1-loop- 2-loop, and 3-loop diagrams. The details of this calculation can be 
found in a; the final result is 



/m(A b ) = -0.21m? 



0.54 + 0.72 log 



A, 



2m j 



9% 
m E 



- 0.70 f-^- 

\m E 



(4) 
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The coefficient Jm in (Q) can be expanded in powers of g by setting g E = 
g 2 T and by substituting the expression (Q) for rn E . The complete free energy 
to order g 5 is then F = (/b + fhi)T\ in agreement with the result obtained 
independently by Kastening and ZhaiH Note that the dependence on the 
arbitrary factorization scale cancels between Je and /m, up to corrections 
that are higher order in g, leaving a logarithm of T/ttie- 

We have calculated the free energy as a perturbation expansion in powers 
of g to order g 5 . We now ask how small a s =^£/(4ir) must be in order for 
the perturbation expansion to be well-behaved Era. If the series is apparently 
convergent, then it can plausibly be used to evaluate the free energy. For 
simplicity, we consider the case rif — 3 and choose the renormalization scale 
fi = 2nT which is the mass of the lightest nonstatic mode. The correction to 
the leading order term of the free energy is a multiplicative factor 1 — 0.9a s + 
3.3a^ 2 + (7.1 + 3.5 loga s )a;? — 20. 8a^ 2 . If the temperature is a few times 
above the transition temperature (~ 200 MeV); e.g., T ~ 350 MeV the series 

has the form 1 - 0.27 + 0.54 + 0.26 - 1.03 H . It is clear that the perturbative 

series is not convergent at this temperature. 

/e and rriE contain the contributions of order T. In the case of m E , 
which is given by the correction to the leading order result form the series 
1 — 0.2Qa s . The next-to-leading-order correction to m% is smaller than the 
leading-order term only if T > 40 MeV. A similar analysis of /e gives T > 67 
MeV. We conclude that the perturbation series for the parameters of EQCD 
are well-behaved provided that T > 70 MeV. 

Expression (JI]) gives /m which is the contribution from the scale T. The 
next-to-leading-order correction is small if A^ is chosen to be approximately 
mE- The next-to-next-to- leading-order correction is smaller than the leading 
order term only if T > 2 GeV. Thus the perturbation series for Jm is well- 
behaved only for temperatures that are much higher than those required for 
the parameters of EQCD to have well-behaved perturbation series. 

This analysis indicates that the slow convergence of the expansion for F 
in powers of -^al can be attributed to the slow convergence of perturbation 
theory at the scale gT. 
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